Thermal T(ls) and Xbi suppression in ^/snn = 2.76 TeV Pb-Pb collisions at the LHC 
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I compute the thermal suppression of the T(ls) and Xbi states in ^snn = 2.76 TeV Pb-Pb 
collisions. Using the suppression of each of these states I estimate the total Raa for the T(ls) state 
as a function of centrality, rapidity, and transverse momentum. I find less suppression of the Xbi 
state than would be traditionally assumed; however, my final results for the total T(ls) suppression 
are in good agreement with recent preliminary CMS data. 



The behavior of nuclear matter at extreme tempera- 
tures is now being studied with the highest collision en- 
ergies ever achieved using the Relativistic Heavy Ion Col- 
lider at Brookhaven National Laboratory and the Large 
Hadron Collider at CERN. For RHIC js]^ = 200 GeV 
Au-Au collisions, initial central temperatures of T ~ 350 
MeV were generated. For current LHC ^snn = 2.76 
TeV collisions one obtains T - 500 - 600 MeV Q] and 
for upcoming full energy runs with ^snn = 5.5 TeV one 
expects To ~ 700 — 800 MeV. At such high temperatures, 
one expects to generate a quark-gluon plasma (QGP) in 
which the formation of quark bound states is suppressed 
in favor of a dcconfincd plasma of quarks and gluons. 

Suppression of quark bound states follows from the fact 
that in the QGP one has Debye screening of color charge 
[2J. Heavy quarkonium has received the most theoreti- 
cal attention, since heavy quark states are dominated by 
short distance physics and can be treated using heavy 
quark effective theory. Based on such effective theories 
of QCD, non-relativistic quarkonium states can be reli- 
ably described. Their binding energies are much smaller 
than the quark mass ttiq ^ Aqcd (Q — c,b), and their 
sizes are much larger than l/mg. At zero temperature, 
since the velocity of the quarks in the bound state is 
small (v <C c), quarkonium can be understood in terms 
of non-relativistic potential models such as the Cornell 
potential which can be derived directly from QCD using 
effective field theory [3]. 

Using such non-relativistic potential models, studies of 
quarkonium spectral functions and meson current corre- 
lators have been performed [4]. The results have been 
compared to first-principles lattice QCD calculations [5] 
which rely on the maximum entropy method [BJ. In re- 
cent years, however, there has been an important theoret- 
ical development, namely the first-principles calculation 
of the thermal widths of heavy quarkonium states which 
emerge from imaginary- valued contributions to the heavy 
quark potential. The first calculation of the leading-order 
perturbative imaginary part of the potential due to glu- 
onic Landau damping was performed by Laine et al. [7J. 
Since then an additional imaginary-valued contribution 
to the potential coming from singlet to octet transitions 
has also been identified [8 . The consequences of such 
imaginary parts on heavy quarkonium spectral functions 



,9J , perturbative thermal widths [7J [S] , and quarkonium 
suppression in a T-matrix approach |10| have recently 
been studied; however, these studies were restricted to 
the case of an isotropic thermal plasma, which is only 
the case if one assumes ideal hydrodynamical evolution. 

The calculation of the heavy quark potential has 
since been extended to the case of a plasma with finite 
momentum-space anisotropy for both the real [11] and 
imaginary [12] parts. Additionally, the impact of the 
imaginary part of the potential on the thermal widths of 
the states in both isotropic and anisotropic plasmas was 
recently studied |13j . The consideration of momentum- 
space anisotropic plasmas is necessary since, for any fi- 
nite shear viscosity, the quark-gluon plasma possesses lo- 
cal momentum-space anisotropies |14Lll5j . Depending on 
the magnitude of the shear viscosity, these momentum- 
space anisotropies can persist for a long time and can be 
quite large, particularly at early times or near the edges 
of the plasma. 

This has motivated the development of a new dynami- 
cal formalism called "anisotropic hydrodynamics" (aHy- 
DRO) which extends traditional viscous hydrodynamical 
treatments to cases in which the local momentum-space 
anisotropy of the plasma can be large [TH [TS] . In this 
Letter I present first results of the combination of the 
aHydro temporal evolution of Ref. [H] together with 
results obtained in Ref. [13j for the real and imaginary 
parts of the binding energy. Using this framework I com- 
pute the suppression of the T(ls) and Xbi mesons as a 
function of centrality, rapidity, and transverse momen- 
tum. 

Model potential: The phase-space distribution of glu- 
ons in the local rest frame is assumed to be given by 

/(X,p) = / iso (VP 2 + £(P ' n ) 2 /Phard) where p h ard IS & 

scale which specifies the typical momentum of the par- 
ticles and can be identified with the temperature in an 
isotropic plasma (£ = 0) [TBJ. In general, the param- 
eter £ measures the degree of anisotropy of the plasma 
via £ = §(pj_)/(P*) — 1 where p z and pj_ are the par- 
tonic longitudinal and transverse momenta in the local 
rest frame, respectively 

The perturbative heavy quark potential as function of 
£ has been evaluated previously and has both real [TT] 
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and imaginary contributions [12] • For N c = 3 and Nt — 
2 the real part of the resulting potential can be well- 
approximated by ^[Vport] = — aexp(— fir)/r with 

(ji_Y 4 _ u J u ^(l + Q 2 (cos(2g)-l) \ 

(i) 

where a = 4a s /3, mfj = (1.4) 2 167ra s Ph ard /3 is the 
isotropic Debye mass, and 9 is the angle with respect to 
the bcamline. The factor of (1.4) 2 accounts for higher- 
order corrections to the isotropic Debye mass |17) . 

However, for describing finite-temperature states 
which can have large radii compared to Aq^ d , one must 
supplement the perturbative short range contribution 
above by a long range contribution. Following previous 
work [TT], I generalize the Karsch-Mehr-Satz potential 
[2] by including the anisotropic mass scale /x given in 
Eq. ([!]) in place of the isotropic Debye mass and adding 
the entropy contribution necessary to obtain the internal 
energy of the states. Such a construction agrees well with 
lattice measurements of the real part of the heavy quark 
potential [TT] . The resulting model for the real part of 
the heavy quark potential is 

mV] = -— (1 + ^r)exp (~nr) + — [1 - exp(-ur)] 
r fi 

s °- 8cr 

-ar exp(-/xr) — , (2) 

m l Q r 

where the last term is a temperature- and spin- 
independent quark mass correction p~8] and a = 0.223 
GeV is the string tension. Here I ignore the effect of the 
running of a s and fix a — 0.385 to match zero tempera- 
ture binding energy data for heavy quark states [TT] . 

For the imaginary part of the model potential I use the 
imaginary part of the perturbative heavy quark potential 
which has been calculated to linear order in £ 

%[V peit ] = -aPhard{0(r)-£[V'i(M)+V>2(M)]}, (3) 

where f = m^r and </>, ipi, and "02 are defined in Ref. |12j . 

The full model potential is given by V = 3?[V] + i%[V] 
and can be used in the Schrodinger equation. To solve 
the resulting Schrodinger equation I use the finite differ- 
ence time domain method [Jj5] extended to the case of a 
complex- valued potential |13j . For the temperature and 
anisotropy dependence of the resulting real and imagi- 
nary parts of the binding energies for the T(ls) and \bi, 
I refer the reader to Ref. [13j . For a point of reference, 
in an isotropic plasma the medium-induced width of the 
T(ls) is approximately SfE'bind] ~ 0.211 T. 
Initial conditions and dynamics: Solution of the 
Schrodinger equation gives the real and imaginary parts 
of the binding energy of the states. The imagi- 
nary part defines the instantaneous width of the state 
9[-E b ind(Phard,£)] = -T T (pbiud, £)/ 2 - However, one must 



account for the complete disassociation of the states when 
^[-Ebind] < 0. I implement this by assigning a large width 
of 10 GeV - (0.02 fm/c)- 1 to states when 5ft[£bmd] < 0. 
The final results are insensitive to the value of this rate, 
as long as it is taken to be greater than 0.5 GeV. Given 
the binding energy data, I evolve the system using the 
non-boost invariant aHydro equations of Ref. [TH] . Us- 
ing the output I can compute the temporal evolution of 
the thermal width which gives IV (r). 

The resulting width IV (t) implicitly depends on the 
initial temperature of the system. I vary the assumed 
plasma shear viscosity to entropy density ratio 4ttt]/S — 
{1,2,3} and for zero impact parameter collisions I use 
central temperatures of T = {520,504,494} MeV, re- 
spectively. The central temperature for Airrj/S = 1 is 
fixed based on the 3+ld viscous hydrodynamical simula- 
tion of Schenke et al which reproduces the particle spec- 
tra and elliptic flow seen in y / i]v r /v = 2.76 TeV Pb-Pb col- 
lisions pQ. The central temperatures at A-n-q/S = {1, 2, 3} 
were chosen in order to keep dN^/dy = 1400 fixed for 
different assumed viscosities. For the rapidity depen- 
dence of the initial temperature, I use the Gaussian ra- 
pidity profile specified in Ref. [15] . Finally, I assume an 
initial anisotropy of £o = 0. Since the current aHy- 
dro implementation does not include transverse dynam- 
ics, I model the transverse evolution at zero and finite 
impact parameter as a set of decoupled longitudinally 
expanding plasmas with initial temperatures given by 
T(xj_,&) = T [AT part (x ± ,6)/iV part (0,0)]3, where the par- 
ticipant density is computed using the Glauber model 
with a Woods-Saxon profile and (Jnn — 62 mb. 

At each transverse point I then evolve the system us- 
ing aHydro starting from tq = 0.3 fm/c, terminating the 
evolution at a final time, rj, when the local energy den- 
sity becomes less than that of an N c = 3 and Nf — 2 ideal 
gas of quark and gluons with a temperature of T = 192 
MeV. At this temperature plasma screening effects are 
assumed to decrease rapidly due to the transition to the 
hadronic phase and the widths of the states will become 
approximately equal to their vacuum widths. 

Formation time: It is also important to consider the 
effect of the time-dilated formation time on the suppres- 
sion of the T(ls) and \bi states [20]. The formation time 
of a state can be estimated by the inverse of its vacuum 
binding energy. Here I take r f ° rm = {0.2,0.4} fm/c for 
the T(ls) and \bi, respectively, in their local rest frame. 
These formation times are consistent with the states' re- 
spective vacuum binding energies. In the lab frame the 
formation time depends on the transverse momentum of 
the state via the gamma factor Tf orm (px) = l T fo X m ~ 
EtT{ otto /M where M is the mass of the relevant state. 
For averaging over transverse momenta I assume that 
both states have a 1/Ej, spectrum which is consistent 
with the high-pT spectrum measured by CDF [21] . 

The suppression factor: Having obtained the 
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FIG. 1. Rapidity- and pr- averaged Raa for T(ls) and \bi as 
a function of N pal t using 4nr//S = 1. 
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FIG. 2. Rapidity- and pr-averaged Raa f° r the T(ls) as a 
function of N ps , Tt for 4nri/S £ {1,2,3}. 
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FIG. 3. Centrality- and pr-averaged -Raa f° r the T(ls) as a 
function of rapidity for 47177/5 £ {1, 2, 3}. 
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FIG. 4. Centrality- and rapidity-averag ed R { Xa for the T(ls) 
as a function of rapidity for 4nri/S € {1, 2, 3}. 



spatial and temporal evolution of the widths I can 
compute the resulting nuclear suppression factor, 
Raa- Starting from the extracted time-, transverse- 
coordinate-, and rapidity-dependent instantaneous de- 
cay rate, I integrate to obtain j(x±,px, <T, b) = 
e(r/ - T im . m {p T ))J^ ax{TtaiirL{pThTo] dTT T (T, X± ,^b), in 
proper time |22j where <; is the spatial rapidity. 
From this I extract Raa via Raa(^-±,Pt-, <Tj b) — 
exp(— 7(x_l,pt, <>, b)). Finally, I average over xj_ tak- 
ing into account the transverse dependence of the hard- 
particle production probability via (Raa(PTi ?> b)) = 
[J x± dx± T A a (x_l ) Raa (xx , Pr, b)}/ [f x± dx± T A a (x_l)] . 

Feed down: Since a certain fraction of T(ls) states 
produced in high energy collisions come from the de- 
cay of excited states, when computing the total Raa 
for the T(ls) one must also consider the suppression of 
the excited states which decay or "feed down" to it. In 
this work I have only computed the Raa for one ex- 
cited state (xbi), so I can only estimate the full feed 
down effect. I will estimate the full feed down effect 
by assuming that all excited states have the same Raa 
as the Xbi (shown in Fig. [I]). With this assumption 



the T(ls) Raa including feed down can be written as 

= a;i?S r ° Und StatC + (1 - X )i?-«ted states where % ig 

the percentage of T(ls) states which are produced di- 
rectly. Measurements of bottomonium state feed down 
in s/s = 1.8 TeV pp collisions at CDF 23J with a cut 
p^, > 8.0 GeV/c find that the percentage of directly pro- 
duced T(ls) states is [50.9 ± 8.2 (stat) ± 9.0 (syst)]%. 
In all plots shown I use x = 0.51. 

RESULTS: In Fig.[T]l show the N paxt dependence of Raa 
for T(ls) and Xbi- As can be seen from this figure, de- 
spite the fact that the initial temperature is not high 
enough to completely dissociate the T(ls), there is still 
a suppression due to the in-medium decay. At these en- 
ergies we see a somewhat similar suppression pattern for 
the Xbi ■ This may seem paradoxical since the naive melt- 
ing temperature for the Xbi is ~ 345 MeV; however, it 
is important to consider the finite formation time of the 
Xbi and the transverse dependence of the temperature. 
In practice, one finds that there are still quite a few Xbi's 
generated from the low temperature regions. After aver- 
aging over the pr-spectrum and geometry, the final result 
is the one shown in the Fig. [T] 
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In Fig. [2]l show the N part dependence of Raa f° r three 
different values of the plasma shear viscosity, r\. The ex- 
perimental data points shown are preliminary data from 
the CMS collaboration [24 . In all cases, the statistical 
error is indicated by the narrow darker (dark red online) 
error bar while the systematic error is indicated by the 
broader (purple online) shaded region. In Fig. [3] I show 
the rapidity dependence of the T(ls) Raa- I n Fig- [f] 1 
show the pt dependence of R 1 Xa f° r 47T7//5 = 1. In this 
figure I averaged over the p-r-bins specified by the exper- 
iment, < p T < 6.5 GeV, 6.5 GeV < p T < 10 GeV, and 
10 GeV < p T < 20 GeV. 

Figures. [2j [3j and [4] taken together demonstrate a 
reasonably good agreement between theory and exper- 
iment; however, the pt dependence of the theoretical 
result seems to be much flatter than the experimental 
results. Practically speaking though, based on the lim- 
ited number of events, it is hard to draw firm conclusions. 
Looking forward it will be important to include the sup- 
pression of the other relevant excited states, e.g. T(2s) 
and T(3s), and the effect of transverse expansion in the 
dynamics. Both of these goals represent work in progress. 
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